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Digital filter and listening device. 

AREA OF THE INVENTION 
5 Traditional digital FIR filters have linear frequency resolution. Thus, these filters are not 
appropriate for matching the frequency resolution of the ear, nor for realizing demands 
for steep filter transitions in the low frequency area. Warped FIR filters can be designed 
to have nonlinear frequency resolution and are thus more appropriate for the described 
tasks. 

10 

In hearing aids and other sound processing devices, where the sound processing 
characteristics changes online, for instance in response to environmental changes 
different types of adaptive filter approaches or compression strategies are typically used 
to change the sound processing characteristics. In adaptive filters, FIR types have the 
15 advantage that they cannot become unstable, due to internal feedback in the filter, as it is 
known from adaptive HR filters. Warped FIR filters do not suffer from the feedback 
problem as long as the warping parameter is numerically smaller than one. 

Designing filters on a warped frequency axis is described by Makhoul and Cosell in 
20 1976 and Strube in 1980. An overview of Frequency- Waiped Signal Processing for 
Audio Applications is given by Harma et al. 2000. 

BACKGROUND OF THE INVENTION 

25 From EP 1 191 814 Al it is known to use a warped filter in a feedback adaptive filter 

approach, with the purpose of reducing howling or feedback, in the hearing aid. Also it is 
known to use the warped filtering technique for other applications, for instance including 
loudspeaker equalization, audio coding and adaptive feedback cancellation. 

30 Further it has been proposed to use warped filtering in an online feed forward filter 
design approach in a hearing aid. 

A regular FIR filter has the transfer function 
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H(z) = ^ ^ (z _I )* delay line elements = z" 7 

By replacing the delay elements with warped delay elements a warped FIR (WFIR) filter 
is obtained 

F (z) = 2* O" 1 )* dela Y line elements = w' 1 
where w= (1 - A*" 1 ) /(z" 1 - X) as shown in Figure 1. 

A warped FIR filter (WFIR) is shown in fig. 1 with warped delay line elements 

w = 

1-Xz 1 

The warping parameter 0<X<\ determines the amount of warping. The filter F(z) can be 
designed using regular FIR filter design techniques, e.g. the Fourier method, except the 
target frequencies are warped. A simple design example is shown in figure 2. Sample 
frequency is 20kHz, X = 0.5. The target filter is a triangular bandpass filter (in absolute 
magnitude) withpassband from 166Hz to 332 Hz. 

The FIR and WFIR filters are designed by similar approaches, the only difference being 
that the WFIR filter is designed on a prewarped frequency axis and the FIR filter being 
designed on a linear frequency axis. The result is easily seen from figure 2. The warped 
filter can match the target better due to more frequency resolution for low frequencies. 
The FIR filter can not match the steep slopes of the target curve at frequencies below 500 
Hz. The filter resolution achieved by the FIR filter at 500 Hz is already achieved in the 
WFIR filter at 167 Hz, when X is set to 0.5 at a 20 kHz sampling frequency. 

It is easily shown, that the frequency resolution of the warped filter is increased by a 
factor determined by the warping parameter, namely the expansion at 0 Hz: 

A/(0) = i±ii 

1 — A 

For X = 0.5 this factor is 3, which means that a warped FIR filter has a resolution at low 
frequencies comparable to a FIR filter which is 3 times longer (at the expense of 
resolution at the high frequencies). 
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The benefit of warped filters is that they can tune their frequency resolution to any 
frequency region needed. If high resolution is needed at high frequencies the X parameter 
must be set at for instance -0.5 relative to 0.5 when high resolution is wanted for low 
frequencies. 

The cost for the increased flexibility is an increase in computational complexity. But 
when selecting X appropriately, for instance to 0.5 the increase in computational 
complexity is low. 

Lower average throughput delay is achieved when matching WFIR and FIR filters of the 
same computational complexity. And the throughput delay is only high at those 
frequencies, where high flexibility is needed. In symmetric FIR filters, the throughput 
delay is constant across frequency. The frequency dependent throughput delay of 
symmetric warped FIR filters is only dependent on the warping parameter X. Thus the 
symmetric WFIR filter has a constant phase, even though the filter coefficients are 
changed (as long as X is kept constant). 

These known approaches do however not take advantage of the full potential of the 
warped filter technique. 

SUMMARY OF THE INVENTION 

According to the invention the warped FIR filter and the FIR filters are combined in an 
attempt to get the "best of both worlds". 

This is achieved with a digital filter comprising at least a first and a second delayed 
summation line, whereby one of the delay lines comprise a warped FIR filter and the 
other line is a FIR filter. According to the invention: 

• at least two delayed summation lines are present, 

• midpoint elements from the first delayed summation line are used as input to the 
second delayed summation line, 
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• output from the first delayed summation line are delayed with a delay 
corresponding to the total delay of the second delayed summation line and 

• the delayed output from the first delayed summation line is added to output from 
me second delayed summation line to form an output. 

The idea is to let the warped filter process one part of the frequency range, for instance 
the lower frequencies, and to let the FIR filter process the remaining frequency range in 
thxs case the high frequencies. This allows the FIR filter to be shorter, and the WFIR 
filter can operate with intensified warping so that it too will be shorter. 

In the z-domain the filter structure according to the invention has the following 
expression: 



y(z) = x(z) 



M-l 



J-l 



L "'=0 J=0 



where the delay lines are tapped halfway 
shown in Figure 4 and Figure 5, i.e.: r = ^* and m- M ~ 1 



as 



Here the summation over w comprises the warping part of the filter and the summation 
over z comprises the FIR part of the filter. 

In an embodiment of the invention the first delayed summation line is a warped 
summation line and the delay of the output of the first delayed summation line is a 
simple delay. 

m another embodiment, the second delayed summation line is a warped summation line 
and the delay of the output of the first delayed summation line consists of warped delay 
elements. 



M yet another embodiment the warped FIR filter comprises a number of warped filter 
sections, whereby midpoint elements of each section are used as input for the next 
section and where further the output of the same sections are delayed for in-phase 
summation with the output from the next section. Using this embodiment allows 
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dedicated filtering of particular parts of the frequency axis defined by the warped filter 
segments. 

The invention also comprises a listening device. The listening device is adapted to 
receive an input signal and has a signal path from the input to a receiver for delivering a 
sound signal to the ear of the user, where at least part of the signal path is digital and 
where the signal path comprises a digital filter as describe above. 

If the listening device is a hearing aid it allows the hearing aid to realize Warped FIR 
filters which have a smaller group delay than similar traditional FIR filters with the samt 
low frequency resolution. The difference is largest for high frequencies, where the delay 
through the digital WFIR filter can be considerably lower than for the FIR filter. In 
hearing aids it is very important to minimize the delay, since delays have a negative 
impact on the sound perception of the user. These negative effects can for instance be 
comb filter effects due to interactions between direct sound through the vent and the 
delayed and amplified sound from the hearing aid. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 Shows a warped FIR filter structure according to the prior art, 

Fig. 2 is a diagram with the transfer functions for warped and regular FIR filters and 

target, 

Fig. 3 shows a simple block diagram of the combined WFIR and FIR solution according 
to the invention, full lines indicate signal flow and dashed lines indicate parameter flow, 
Fig. 4 is an example of parallel FIR (of order 2k) and WFIR filter (of order 2p) with 
matched phase, 

Fig. 5 is an example of parallel FIR (of order 2k) and WFIR filter (of order 2p) with 
matched phase, advantageous for special cases, for instance for small warped filter 
lengths, 

Fig. 6 shows a combination of 3 WFIR blocks and 1 FIR block 
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The use of the combined FIR and WFIR filter according to the invention in a feed 
forward adaptive filter approach is shown in figure 3. The controller feeds filter 
coefficients to the Warped FIR filter and the FIR filter online, and thus changes the 
characteristics of the hearing aid signal processing based on at least changes in the 
digitized input signal. 

The filter structure in fig. 4 is designed according to the invention, z" 1 indicates a delay 
element with a delay of one sample, z~ k indicates a delay element with a delay of k 
samples, and w~ x indicates a warped delay element, with ^(l-^- 1 ) _ x) 
The warped delay line is tapped mid-way and fed as input to the FIR filter and a simple 
delay, matching that of the FIR filter, is placed at the output of the WFIR filter. This 
causes the two signal paths to match each other in phase, so that the adder at the output 
will not cause signals travelling the two paths to cancel out. 

In some cases (for small warped filter lengths), it may be an advantage to let the FIR and 
warped FIR filters trade places. This is illustrated in Figure 5 where z~ ' indicates a 
delay element with a delay of one sample, w~ k indicates a warped delay line with k 
warped delay elements, andM*=(l- Az~ l ) /(z~ l - X) . 

In stead of having only one warped filter section, it is also possible to extend the 
combined WFIR and FIR structure, so that multiple WFIR sections with different 
warping parameters X is used to address different frequency regions. An example of such 
a solution is shown in figure 6. 

The figure shows a combination of 3 WFIR blocks and 1 FIR block. This combination 
allows the use of 3 different warping parameters X, which allows the WFIR filters to 
have different frequency resolution, allowing to match very specific needs on 3 different 
warped frequency scales, as well as on the traditional linear frequency scale. 



